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CALIBRATION  OF  THE  OPTICAL  QUALITY  OF  A  LARGE  TELESCOPE  SYSTEM 
Wu  Yi , Wang  Yingj ian, Wu  Xubin 


ABSTRACT  Use  is  made  of  optical  correlation  functions 
associated  with  atmospheric  trubulence  to  precisely  determine 
turbulence  optical  correlation  times.  In  conjunction  with  this 
within  these  time  periods — after  atmospheric  turbulence  effects  are 
corrected  through  imagery  system  optical  transfer  functions,  large 
telescope  (500mm  aperture,  500m  focal  length)  system  optical 
quality  is  calibrated. 

KEY  WORDS  Correlation  function  Turbulence  Large  telescope 
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0  INTRODUCTION 

In  experimental  research  during  actual  atmospheric  laser 
transmission  or  imaging  of  celestial  bodies  through  the  atmosphere, 
it  is  often  necessary  for  relatively  large  aperture  and  long  focal 
length  telescopes  to  be  transmission  or  reception  systems.  The 
purposes  are:  1)  to  increase  beam  transmission  efficiency  in  order 
to  reduce  beam  transmission  diffraction  angles;  2)  to  increase 
spacial  resolutions  associated  with  imaging  systems.  However,  not 
only  atmospheric  turbulence  effects  will  restrict  their 
realization.  Optical  gualities  of  telescopes,  themselves  .are  also 
an  influencing  factor.  The  reason  is  that  it  is  very  difficult  for 
large  telescope  systems  to  reach  ideal  levels  that  is,  diffraction 
limitation  systems.  As  a  result,  during  the  carrying  out  of 
experimental  research  using  actual  trubulent  atmosphere  as  a 
transmission  medium,  there  is  a  need  to  carry  out  calibrations  on 
telescope  system  optical  qualities.  Normally,  ideal  system 
diffraction  limitation  multiples  are  used  as  characteristic. 
However,  speaking  in  regard  to  large  telescope  systems— in 
particular,  long  focal  length  systems— it  is  very  difficult,  within 
the  limited  space  of  the  laboratory,  to  carry  out  optical  quality 
calibrations.  There  is  a  need  to  carry  out  this  work  in  the  actual 
atmosphere.  Using  the  real  atmosphere  as  background— when  carrying 
out  calibration  of  large  system  optical  qualities— it  is  then 
necessary  to  experience  influences  associated  with  turbulence 
effects.  As  a  result,  there  is  a  need  to  carry  out  corrections 
with  regard  to  turbulence  effects.  Actual  atmospheric  turbulence 
often  shows  itself  as  a  nonsteady  random  process.  However,  with 
respect  to  descriptions  associated  with  turbulence  effect  optica 
parameters,  their  source  is  then  steady  random  processes  [1/2]. 
This  article  makes  use  of  correlation  functions  associated  with 
finite  sequence  random  sampling  in  regard  to  turbulence  effect 
optical  parameters.  First  of  all,  precise  determinations  are  made 
of  correlation  times  associated  with  atmospheric  turbulence .  After 
that  within  turbulence  correlation  time  periods,  long  exposure 
statistics  [1,2]  are  made  with  respect  to  atmospheric  turbulence 
optical  transfer  functions  «OTF..  Finally,  calibrations  are 
achieved  of  large  telescope  system  diffraction  limitation  multiples 
after  turbulence  effect  corrections.  The  actual  research  object  in 
this  article  is  a  500mm  aperture  variable  focal  length  Cassegrain 
telescope  system  developed  by  the  Chinese  Academy  of  Sciences 
Chengdu  Photoelectric  Institute. 

1  EXPERIMENTAL  PRINCIPLES 

Fig. i  is  a  telescope  system  calibration  experiment  light  path. 

The  telescope  has  focus  adjusted  to  500m.  CCD  is  used  in  recording 
light  source  imagery  point  spread  functions  and  beam  arrival  angles 
calculated  under  the  influence  of  turbulence  effects.  /511 
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Fig.l  Experiment  Optics 
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1.1  Telescope  Diffraction  Limitation  Multiples 

Taking  telescopes  and  atmospheric  turbulence  together  as  an 
integrated  imaging  system  forming  images  with  regard  to  laser  light 
sources  at  places  500  meters  distant  (approximating  a  point  light 
source),  light  striation  point  spread  functions  on  CCD  image 
surfaces  at  this  time  are: 

PSF  =  PSFTur*PSFTel  (1) 

In  this,  PSFTur  is  the  atmospheric  turbulence  point  spread 

PSFTel  is  the  telescope  system  point  spread  function.  *  stands  for 
scroll  product.  The  modulation  transfer  function  is  the  Fourier 
transform  mode  [3]  associated  with  point  spread  functions,  that  is 

MTF  =  |OTF  |  =  F(PSF)  (2) 


F(  )  stands  for  Fourier  transform.  Therefore 


MTF 


Tel- 


MTF 

MTFTttr 


(3) 


MTF  is  obtained  by  the  Fourier  transform  of  imagery  surface  point 
spread  functions.  During  long  exposure  statistics  by  turbulence 
optical  statistics  theory,  MTFTur  is  then  [2] 

MTFXur  =  exp[— 3.44(  —  )5/3u3/3]  (4) 


In  this,  D  is  telescope  diameter,  ro  is  atmospheric  turbulence 
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coherence  length.  u=fkA/D  is  normalized  spacial  frequency,  f  is 
focal  length,  k  is  spacial  frequency.  A=0.6328Mm.  Ideal 
telescope  imaging  system  modulation  transfer  functions  are  [2] 

MTFt,=  -j  (cos-'u-uVl-u2  )  (5) 

The  cut  off  spacial  frequency  is  kc(u=l )=D/(f X ) .  Resolving  power 
is 

r  ^ 

Rq-4(-j  )2j  «(cos-'«-«Vl-u2  )d“  (6) 


Actual  telescope  system  spacial  resolution  is  [2] 


uMTFTel(u)du 


(7) 


Therefore, 

precisely 


telescope  system  diffractive  limitation  multiples  can  be 
determined  from  the  equation  below 


B 


2j  u(cos ~‘u  — u-J  I  —  u2  )du 


-f 


«MTFTeI(«)du 


} 


1/2 


(8) 


Jo 


1.2  Precise  Determination  of  Atmospheric  Turbulence  Correlation 
Time  Periods  /512 

We  know  that  theoretical  descriptions  with  regard  to 
atmospheric  turbulence  optical  parameters  are  based  on  steady 
random  processes  [1] .  Due  to  the  fact  that  actual  atmospheric 
turbulence  is  often  a  nonsteady  random  process,  it  is  necessary,  as 
a  result,  to  carry  out  statistical  descriptions  with  respect  to 
atmospheric  turbulence  optical  parameters  under  a  presupposition 
finite  steady  random  processes.  It  is  generally  believed  that  this 
operation  should  be  carried  out  within  turbulence  effect 
correlation  time  periods.  Fluctuations  in  beam  arrival  angles  and 
their  correlation  functions  have  the  relationship  below  [1] 


B(  x)e-io'dx= 


d  t 


1/2 


(9) 


In  this,  B(r)  is  the  turbulence  optical  correlation  function. 
fa(t)  is  a  time  function  associated  with  oscillations  in  beam 
arrival  angles.  The  equation  above  clearly  shows  that  power 
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spectra  associated  with  oscillations  in  arrival  angles  are  equal  to 
Fourier  transforms  of  correlation  functions.  This  article  makes 
use  of  sampling  sequences  associated  with  light  striations  on  CCD 
imaqe  surfaces .  Calculations  are  made  of  beam  arrival  angle  time 
period  sequences/  and,  in  conjunction  with  that,  discrete  Fourier 
transforms  are  carried  out  on  them.  After  that,  use  is  made  of 
discrete  forms  of  equations  above  to  carry  out  Fourier  inverse 
transforms,  obtaining  the  atmospheric  turbulence  optical 
correlation  function  B(r).  It  is  generally  believed  that,  when 
B(t)  drops  down  to  be  zero,  the  corresponding  time  to  it  the 
turbulence  correlation  time. 

2  EXPERIMENTAL  PARAMETERS  AND  RESULTS 


Fig. 2 (a)  Power  spectrum  of  arrival-angles 


Fig. 2(b)  Correlation  function  from  power  spectrum  of  arrival 


u 


Fig. 3  MTFs  with  normalized  spatial  frequence 


In  accordance  with  Fig.l,  a  1mm  He-Ne  small  aperture  laser  is 
shot  out  through  500m  of  turbulent  atmosphere  and  received  at 
telescopes.  Through  lenses,  light  striations  form  imagery  on  CCD 
image  surfaces.  At  this  time,  due  to  the  fact  that  light  source 
small  aperture  dimensions  are  already  smaller  than  500mm  aperture 
ideal  systems  associated  with  limit  resolution  scales  at  500m. 
distances,  it  is,  therefore,  possible  to  approximate  their  being 
point  light  sources.  Optical  striation  point  spread  functions 
associated  with  CCD  image  surfaces  are  desribed  by  equation  (1). 
Lens  focal  length  is  450mm.  During  experiments,  at  the  same  time, 
option  is  made  for  the  use  of  two  CCD  units  to  record  optical 
striation  imagery.  The  sampling  frequency  is  380Hz.  Continuous 
recording  of  8192  frames  is  used  in  obtaining  arrival  angle 
oscillation  power  spectra  and  atmospheric  turbulence  optical 
correlation  functions.  Due  to  the  fact  that  the  sampling  frequency 
is  fast,  it  is  possible— after  going  through  high  speed  Fourier 
transforms  of  8192  discrete  data — to  obtain  arrival  angle 
oscillations  power  spectra  associated  with  frequencies  that  are 
0.045-190Hz.  At  the  same  time,  from  variances  associated  with 
arrival  angle  oscillations,  it  is  possible  to  obtain  turbulence 
lateral  coherence  length  [3,4] 

r0=  0.545D-,/5l6/s</«2>  3/5  (10) 


fa  is  arrival  angle.  Experimental  value  associated  with  ro  is 
90mm  Fig.  2(a)  is  power  spectra  associated  with  arrival  angle 
oscillations.  Fig.  2(b)  is  atmospheric  turbulence  normalized 
correlation  functions  solved  for  from  discrete  forms  associated . 
with  equation  (9)  on  the  basis  of  Fig.  2(a)  data.  From  the  Fig.  s 
in  question,  it  is  possible  to  see  that  at  11  seconds  correlation 
function  curves  show  the  appearance  of  negative  values  (not  due  to 
measurement  fluctuations )— ~that  is,,  negative  correlation. 
Therefore,  turbulence  correlation  time  is  selected  as  11  seconds. 
The  sampling  frame  frequency  associated  with  the  second  CCD  unit  is 
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10Hz.  However,  the  image  element  resolution  is  higher  than  the 
former.  It  is  used  in  obtaining  imagery  system  modulation  transfer 
functions  MTF.  From  a  summation  of  110  frames  of  optical  striation 
imagery  within  11  seconds,  average  FFT  obtain  long  exposure  MTF. 
From  equation  (4),  calculations  are  done  of  atmospheric  turbulence 
modulation  transfer  functions  MTFTur  in  long  exposure  situations. 
Fig. 3  gives  curves  associated  with  long  exposure  MTF  and  long 
exposure  MTFTur.  Finally,  from  equation  (3),  solutions  are  made 
for 

telescope  system  MTFTel — after  long  exposure  average  turbulence 
effect  corrections — as  shown  in  Fig. 3.  From  equation  (8),  the 
telescope  system  diffraction  limitation  multiple  B=3.087  is 
obtained. 

3  CONCLUSIONS 

In  cases  where  corrections  of  the  influences  of  atmospheric 
trubulence  effects  are  considered,  making  use  of  square  root 
reciprocal  relationships  associated  with  resolution  specific  values 
between  optical  imagery  system  spacial  resolution  and  ideal  system 
spacial  resolution  to  precisely  determine  large  telescope  system 
diffraction  limitation  multiples  is  one  effective  means.  In  one 
area,  it  assists  in  our  understanding  optical  qualities  of  actual 
telescope  systems  in  order  to  quantify  performance  indices 
associated  with  sending  systems  as  supplied  by  laser  atmospheric 
transmission  tests.  In  another  area,  the  difficulties  associated 
wj_th  carrying  out  indoor  calibration  of  large  telescope  systems 
were  resolved,  and,  in  conjunction  with  that,  methods  were  supplied 
to  eliminate  atmospheric  turbulence  effects. 
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MEASUREMENT  OF  WATER  VAPOR  CONTENT  IN 
SINGLET  OXYGEN  GENERATOR  BY 
SPECTROMETRIC  METHODS 

ABSTRACT  Water  vapor  content  is  one  of  the  important 
performance  parameters  associated  with  chemical  singlet  oxygen 
generators  (SOG).  Because  it  not  only  has  a  quenching  effect  on 
02  (1A)  but  also  has  a  severe  quenching  effect  with  regard  to  I* 
in  oxygen  iodine  chemical  laser  (COIL)  optical  cavities,  as  a 
result,  the  measurement  and  control  of  water  vapor  content  at  SOG 
exit  locations  has  important  significance.  The  experiments  in 
question  opt  for  the  use  of  optical  spectrum  strength  comparison 
methods  to  measure  water  vapor  content  amounts  in  SOG.  Making 
use  of  two  dimensional  scanning  and  accumulating  functions  of 
OMA-3,  it  is  possible  to  carry  out  simultaneous  scanning  and 
accumulation  with  respect  to  optical  spectra  at  different 
measurement  points.  Thus,  it  is  possible  to  carry  out  water 
vapor  amount  determinations  with  regard  to  a  number  of  individual 
SOG  measurement  points  at  the  same  time.  Test  measurement 
precisions  can  reach  ±10%. 


KEY  WORDS  Water  vapor  Singlet  oxygen  generator 
Spectrometric  method 

Cui  Tieji  Zhang  Yunlu  Min  Xiangde  Sang  Fengting  Sun  Yizhu 
Chen  Wenwu  Pang  Jingke  Jin  Yuqi  Yang  Bailing 


0  INTRODUCTION 


The  overall  level  associated  with  oxygen  iodine  lasers 
(COIL) [ 1 ]  shows  up  primarily  in  performance  levels  associated 
with  singlet  oxygen  generators  (SOG)  and  media  energy  sources  of 
the  systems  in  question,  adopted  efficiencies  associated  with 
laser  cavity  structure  with  regard  to  laser  media  energies,  and 
12  generator  performance.  As  far  as  02 (1A)  concentrations 
produced  in  chemical  SOG,  which  COIL's  opt  for  the  use  of,  are  - 
concerned,  water  vapor  and  residual  oxygen  content  are  primary 
indices  measuring  SOG  performance.  This  article  will  make 
additional  discussions  with  respect  to  water  vapor  amount 

measurement  methods.  .  ..  , 

Although  there  are  a  relatively  large  number  of  methods 
to  measure  water  vapor  content,  most  of  them  are  comparatively 
difficult.  Primarily,  they  include  absorption  spectrometric 
methods — for  example,  spectrophotometric  methods  [3],  and  diode 
laser  spectrometer  methods  [4].  The  former  makes  use  of  20m  long 
absorption  pools  with  wave  lengths  of  5900nm  and  time  period 
response  differences.  Moreover,  the  instrumentation  equiment 
associated  with  the  latter  is  very  expensive.  However,  its  test 
measurement  precisions  are  relatively  high,  and  it  is  capable  of 
acting  as  contrast  standard.  There  is  also  one  type  that  is  a 
hygroscopic  solid  hygrometer  method  [5].  It  requires  that 
poisonous  gases  cannot  be  contained  in  gas  flows.  As  a  result, 
it  is  very  difficult  to  apply  it  to  COIL  in  SOG.  Comparatively 
speaking,  the  best  method  is  to  make  use  of  comaparative  /516 
methodsassociated  with  spectral  line  strengths  [4].  It  is 
relativelysimple.  Response  is  fast,  and  it  is  comparatively  easy 
to  make  use  of  on  site.  This  article  opts  for  the  use  of  optical 
multichannel  analyzers  OMA.  They  are  capable  of  simultaneously 
measuring  spectral  line  strengths . 
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Fig.l  Schematic  of  cw  supersonic  COIL  setup 


These  experiments  opt  for  the  use  of  experimental  apparatus 
like  that  shown  in  Fig.l.  In  rotating  SOG,  first  of  all, 
premixed  basic  hydrogen  peroxide  BHP  solution  is  added  in.  This 
can  be  prepared  in  accordance  with  fixed  proportions  from  aqueous 
solutions  of  metallic  hydroxides  MOH  such  as  Li,  Na,  K,  and  so 
on,  and  H202  (concentration  generally  is  35-90%).  In  this  way — 

in  BHP  solutions - there  necessarily  exists  quite  an  amount  of 

liquid  water.  In  conjunction  with  this,  when  C12  is  added  in,  on 
BHP  liquid  film  surfaces  which  are  attached  onto  rotating  disks 
and  continually  renewed,  intense  exothermic  chemical  reactions 
are  generated,  producing  02 (‘A)  .  At  the  same  time,  H20  is  also 
produced 

:  H202  +  2MOH  +  Cl2  —  02('A)  +  2MC1  +  2H20  ( 1 ) 


Large  amounts  of  reaction  heat  will  cause  precooled  BHP  to 
violently  rise  in  temperature,  thereby  producing  a  certain  amount 
of  H20  and  H202  vapors  as  well  as  small  liquid  droplets  and 
dilute  gases  flowing  out  together.  Although,  when  droplets  and 
water  vapor  go  through  SOG  exit  cold  traps,  the  majority  of  them 
can  be  cooled  out,  from  the  time  of  their  formation  onward, 
however,  there  will  then  be  certain  gas  phase  and  liquid  phase 
quenching  effects  with  regard  to  02 (lA)  and  so  on. 

In  COIL  activation  media,  water  vapor  is  also  very 


harmful[2].  It  is  not  only  capable  of  making  the  speed  of  12 
molecule  dissociation  drop.  It  is  also  capable  of  making  I* 
produce  severe  gas  phase  quenching.  Going  a  step  further,  COIL 
performance  is  made  to  drop — even  to  the  point  of  being  unable  to 
produce  light.  For  this  reason,  it  is  necessary  to  strictly 
measure  and  control  water  vapor  amounts  in  SOG. 


1  MEASUREMENT  PRINCIPLES  AND  EXPERIMENTS 


In  gas  flows  at  SOG  exits,  concnetrations  of  02(1  A  )  are 
capable  of  reaching  50-90%.  In  gas  flows — at  the  same  time — 
there  also  exist  small  amounts  of  H20,  H202,  C12,  02 (IX),  as  well 
as  02(31),  and  other  such  gases.  During  the  transmission 
process,  they  will  give  rise  to  the  reactions  set  out  below 


02(‘A)  +  02(‘A)  -  02('Z)  +  02(3Z) 


(2) 


K1  =  2.0x10-17  cm3  /  molecule  •  s,  Reference  [6], 


02('A)  +  02(lA) -*■  202(3X)  +  /iv  ,  634nm 


(3) 


K2  =  4.4x10-23  cm3  /  molecule  •  s.  Reference  [7], 


02('I)  +  H20  —  02('A)  +  H20 


K3  =  5.5x10-12  cm3  /  molecule  v6,34  s,  Reference  [6]. 


(4) 
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02('S)  — 02(3Z)  +  bv  ,  762nm 


(5) 


A1  =  0.077s-!  ,  Reference  [8]. 


Reaction  (2)  is  a  formula  for  the  reaction  producing  02 (IS). 


The  production  rate  is 


»,  =  £,[ 02('A)]2 


(6) 


Reaction  (4)  represents  the  quenching  path  associated  with 
•  qnr  in  reality,  there  are  such  collision  quenching 

02 (IS)  m  s  *  H202  C12  02 (3S),  He,  and  so  on.  The  ratios 

companions  as  H20,  H202,  u/,  >  _nect  to  o2(lS)  are 

of  their  cruenching  speed  constants  [9]  with  respect  to  u^li  j 
7x105^3x104^200 : 4 : 1 .  When  temperatures  are  approximately  27 3K, 

SSBSFK^ 


03  =  K3[O2('S)][H2O] 


(7) 


Going  a  step  further, 
dynamic  equilibrium  form  as 


it  is  possible  to  write  out  the 
sociated  with  02 (IS) 


d[02('S)]/dr  -  K,[02('A)]2-  X3[02('S)][H20]  -  Ax[  02('S)] 

(8) 


sar-a^iSt  srs.s.TST'  • . 

and  A1  -  0 . 077S-1  .  One  then  has  >>  A1. ^Therefore, ^ 

^n^Z^tlrs^r^h'^nfu^,  then,  d[02(l„  ]/**>. 

After  rearrangement,  one  has 


|3 


[H10]  =  K,[0!(i4)],/K![01(,2)1 


(9) 


Luminescence  strengths  measured  by  OMA— 3  probes  for  the  two 
reactions  (3)  and  (5)  respectively  are 

/634=^4«2[02('A)]2  (10) 

hu~  (11) 


B  is  the  detection  system  spectral  reception  efficiency  and 
sensitivity  parameter.  Use  is  required  of  standard  light  sources 
or  other  mthods  for  standardization.  Taking  (10)  and  (11)  and 
substituting  into  (9),  one  has 

[H;0]  =  {Kl/K2}{Ai/ K3}  { Pj6il  Pm}  {hulheii  ( 12 ) 


Letting  m={Kl/K2}{Al/K3}{p7 62/p 634},  then 

[H20]  =  m{/634//762}  ( 13 ) 

m  is  overall  system  standardization  constant.  It  is  possible  to 
obtain  this  from  relatively  precise  standardized  methods  of  doing 
test  measurements  on  amounts  of  water  vapor  contained. 

Therefore,  this  method  is  simple  and  of  practical  use.  However, 
its  physical  significance  is  not  as  clear  as  that  of  equation 

( 12 )  • 

If  use  is  made  of  (2),  (4),  and  (5),  as  well  as  the  reaction 
02 (1A)  ->  02 (3S)  +  hv  ( 1270nm) — going  through  the  same  type  of 
handling  processes  described  above — it  is  also  possible  to  obtain 
similar  results  to  equation  (12) 

[H2O]cc(/I270)2/J7tt  (14) 

In  Fig.l,  the  mutual  relationships  between  the  measurement 
system  in  question  and  COIL  are  clearly  shown.  Optical  radiation 
associated  with  reactions  (3)  and  (5)  collected  from  test 
measurement  apertures  is  sent  through  optical  fibers  into 
spectral  instruments  associated  with  optical  multichannel 
analyzer  OMA-3  systems  and  shined  onto  narrow  slits.  In  two 
dimensional  probes  OMA-3  is  capable  of  carrying  out  simultaneous 
scanning  and  accumulation  with  respect  to  two  spectral  lines, 
thereby  improving  overall  measurement  time  response  and  signal  to 
noise  ratios.  Due  to  limiations  associated  with  OMA-3  probe 
spectrum  response  ranges,  it  is  only  possible  to  make  use  of 
equation  (12)  to  measure  water  concentrations.  Fig. 2  is  scanning 
spectra  for  different  amounts  of  water  content  associated  with 
certain  instants  when  C12  is  poured  into  SOG.  Using  the 


intearation  surface  area  ratio  S634/S762  for  wave  lengths 
associated  with  the  two  spectral  lines  634nm  and  762nm  to  act  as 
optical  strenth  ratio  in  equation  (12),  equation  (13)  then 
becomes 

[H20]  =  mSfsJSw  15) 
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(a)Lower  water  content 


(b)Higher  water  content 


Fig. 2  The  spectrum  under  different  water  content 

Under  normal  experimental  conditions,  02 (1A)  concentration 
changes  produced  by  SOG  in  periods  before  and  aft®r  ®xperimen 
are  not  large.  The  lower  [H20]  is,  the  larger  spectral  line 
surface  areas  associated  with  762nm  are — as  is  shown  in  Fig. 2(a). 
wSen  [H2^  is  very  low,  the  leading  role  of  H20  among  collision 
quenching  companions  of  02 (IS)  is  also  gradually  yielded  up  to 
other  companions.  That  is  nothing  else  than  to  say  that  it  is 
only  necessary  that  KM[M]  «  K3[H20]  ,  that  is,  when 

[H,0]»  K„[M]/K,  (16) 

and  then— and  only  then— are  the  conditions  discussed  above  for 
ianorinq  those  quenching  companions  satisfied.  For  example,  when 
in  2he  273K  rtnge,  the  saturation  pressure  of  the  second  main 
quenching  companion  H202  is  0.3  Torr.  Selecting  appropriate 
ouenchinq  speed  constants  [9],  from  equation  (16),  it  is 
^  'Kin  fn  water  vapor  partial  pressures  much  greater 

rha^o  015  Torr  If  calculations  are  done  in  accordance  with  10 
S!d  water  partial  pressures  should  be  at  least  >0.15  Torr. 

Only  at  this  time  are  the  results  obtained  by  equation  (12) 

comparatively  Precise  increases„as  ie  shoWn  iD  Fig.2(b)- 

then  splrtra!  line  areas  are  gradually  smaller  until  they 

arJ  subverted  in  the  signal  to  noise  ratio.  At  this  time, 

measurement  ^rors  enlarge  c^Unls  Silftten  be 

SCrert??e!  ^hi’s  r“soS,  it  is  possible  to  obtain 
suitable  ranges  for  the  test  measurement  methods  in  question 


15 


which  should  be  in  the  interval  where  water  partial  pressures  are 
0.1-5  Torr. 

2  RESULTS  AND  DISCUSSION 

SOG  developed  in  previous  periods  opted  for  the  use  of  tube 
type  dry  ice  cold  traps.  After  adequate  BHP  solutions  were 
added,  it  was  possible  to  make  3-4  iterations  of  tests.  Each 
iteration  was  2-4  seconds.  The  iterval  between  iterations  was 
around  1-10  minutes.  Average  water  content  results  for  three 
continuous  iterations  of  tests  were,  respectively,  0.88,  1.05, 
and  1.30  Torr.  Later,  SOG  structures  were  improved.  Option  was 
made  for  refrigerator-alcohol  cold  traps.  In  conjunction  with 
this,  working  under  certain  types  of  operational  methods,  average 
water  content  amounts  associated  with  three  continuous  iterations 
of  tests  were  0.36,  0.24,  and  0.19  Torr.  With  liquid  nitrogen 
(LN2)  cold  traps,  option  is  made  for  certain  methods  of  pouring 
LN2.  Results  were  0.35,  0.60,  and  0.74  Torr.  From  among 
experimental  results,  it  was  discovered  that  LN2  cold  traps  have 
better  cooling  results  than  dry  ice  cold  traps.  Moreover,  water 
content  amounts  associated  with  test  measurements  for  LN2  and  dry 
ice  cold  traps  follow  along  with  the  number  of  experimental 
iterations  and  drop.  This  is  because  cold  trap  efficiencies  are 
in  the  process  of  falling.  With  refrigerator-alcohol  cold 
traps,  water  content  amounts  follow  along  with  the  number  of  test 
iterations  and  go  down.  This  is  because  cold  trap  efficiencies 
are  in  the  process  of  going  up.  It  was  also  discovered  that 
there  is  an  intimate  relationship  between  such  things  as  water 
content  amounts  associated  with  the  latter  periods  of  the  same 
experimental  iteration  ,  average  water  content  amounts  and  SOG 
structures,  flow  amounts  and  operating  conditions,  etc. 

The  test  measurement  methods  in  question  make  use  of  optical 
fibers,  taking  optical  radiation  associated  with  measured  points 
and — using  the  same  type  of  conditions — introducing  OMA-3  probes . 
To  this  end,  measurement  instruments  can  be  quite  far  away  from 
COIL  equipment  measurement  points.  Because  of  OMA-3,  it  is 
possible  to  carry  out  simultaneous  scanning,  accumulation,  and 
sampling  of  different  spectra  at  the  same  location,  thereby 
avoiding  using  two  sensors  and  errors  which  are  brought  in  when 
two  light  paths  are  measured.  On  the  basis  of  OMA-3  two 
dimensional  scanning  functions,  it  is  possible  to  take  optical 
radiation  associated  with  the  same  location  (for  example,  in 
front  of  and  behind  cold  traps  or  other  important  locations  to 
study),  and,  using  optical  fiber,  introduce  OMA-3  probes  at  /519 
ifferent  sites  at  the  same  time.  Scanning  and  accumulation  is 
done  with  regard  to  optical  radiation  signals  associated  with 
different  sites  (time  intervals  are  capable  of  reaching  a  few 
tens  of  milliseconds).  Thus,  it  is  possible  to  measure  water 
vapor  amounts  contained  at  various  measurement  points  at  the  same 
time.  Going  a  step  further,  one  obtains,  for  instance,  such 
results  as  cold  trap  efficiencies  and  so  on,  supplying  data  for 
evaluations  with  regard  to  SOG  performance  as  well  as  COIL 


lsvcls  •  • 

From  the  principles  of  test  measurement  methods,  it  is 

nnccihle  to  know  that  spectrometric  methods  are  not  able  to 
measure  water  in  condensed  states  which  is  capable  of  existing  in 

SOG. 

Acknowledgements:  Participating  in  this  experimental  research, 

there  were  also  such  comrades  as  Shen  Huihua,  Li  Mingsheng,  Wang 
Dongmei,  Liao  Weiguang,  Cong  Jiagang,  and  so  on.  For  this,  we 

express  our  thanks. 


REFERENCES 


1  ft***.  1993,  5(3):389 

2  Zagidullin  M  V  el  al.  Sov  J  Quantum  Electron,  1987,  17(3):320 

3  Bonnet  J  et  al.  Appl  Phys  Lett,  1984,  45(10):1009 

4  Barnault  B  et  al.  Eighth  International  Symposium  on  Gas  Flow  and  Chemical  Lasers, 
SPIE,  1990,  1397:231 

5  Georges  E  et  al.  ONER  A,  BP72, 92322  Chatillon  Cedex,  France,  1994,  Private  communication 

6  Heidner  R  F  IH  et  al.  J  Phys  Chem,  1983,  87(13):2348 

7  Fisk  G  A  et  al.  J  Chem  Phys,  1983,  77(10):4965 

8  Heidner  R  FI  el  si,  /  Chem  Phys,  1981,  74(10):5618 

9  dtSC:  1992,183 


EFFECT  OF  ATMOSPHERIC  TURBULENCE  ON  PRECISION  OF 

LASER  TRACKING  SYSTEM 
Fan  Chengyu  Song  Zhengfang 

ABSTRACT  This  article  discusses  the  influences  of 
atmospheric  turbulence  on  laser  tracking  system  angular 
precision.  Use  is  made  of  phase  screen  methods  to  calculate 
errors  given  rise  to  by  atmospheric  turbulence  with  regard  to 
laser  tracking  systems.  Calculations  clearly  show  that,  under 
normal  conditions,  tracking  errors  given  rise  to  by  atmospheric 
turbulence  are  capable  of  reaching  several  tens  of  microradians. 


KEY  WORDS  Intensity  fluctuation 
Tracking  error 


Phase  fluctuation 


0  INTRODUCTION 


In  photoelectric  systems  where  there  is  a  need  to  precisely 
determine  such  things  as  target  position,  long  range  laser 
collimation  associated  with  angular  coordinates  [1],  as  well  as 
tracking  systems,  and  so  on,  use  is  normally  made  of  laser 
tracking  detection  systems.  In  References  [3]  and  [4],  the 
tracking  errors  have  already  been  discussed  in  detail.  However, 
only  tracking  errors  given  rise  to  by  dark  electric  current, 
electronic  circuit,  thermal  noise,  and  phase  fluctuations  were, 
principally  involved.  It  was  believed  that  the  influences  given 
rise  to  by  intensity  fluctuations  created  by  atmospheric 
turbulence  were  certainly  not  large.  However,  during  actual 
applications,  it  was  discovered  that  the  influences  of  intensity 
fluctuations  were  very  severe.  In  particuler,  in  situations 
associated  with  low  angle  of  elevation  real  time  tracking, 
tracking  systems  would  often  miss  targets.  Investigating  the- 
causes  carefully,  they  are  possibly  due  to  the  fact  that,  when 
solving  for  angular  target  miss  amounts,  the  intensity 
distributions  associated  with  optical  faculae  themselves  are 
assumed  to  be  uniform  or  presenting  a  Gaussian  distribution.  In 
the  four  quadrants,  differences  in  light  intensities  are  only  due 
to  distribution  shares,  that  is,  angular  positions,  being 
different.  However,  in  reality — besides  giving  rise  to  phase 
fluctuations — atmospheric  turbulence  will  also  create  light 
intensities  making  random  distributions.  This  article  will 
comprehensively  consider  the  influences  of  phase  fluctuations  and 
light  intensity  fluctuations  with  regard  to  laser  tracking  system 
angular  measurement  precisions  when  beams  are  transmited  through 
trubulent  atmosphere.  Initial  results  clearly  show  that  the 
influences  of  turbulent  atmosphere  with  respect  to  tracking 
errors  will  be  larger  than  noise  associated  with  detectors 
themselves.  In  particular,  the  influences  of  optical  intensity 
fluctuations  are  larger.  Under  normal  circumstances,  tracking 
errors  are  capable  of  reaching  several  tens  of  microradians. 


1  DETECTION  PRINCIPLES 


Option  is  made  for  the  use  of  four  quadrant  photoelectric 
diode  laser  tracking  principles  as  shown  in  Fig.l.  Assuming  that 
Gaussian  beams  are  transmited  through  atmospheric  turbulence  to 
photoelectric  detection  devices  as  shown  in  Fig.l,  if  the 
coordinate  origin  point  is  selected  at  the  center  of  incident 
light  faculae,  incident  optical  intensity  distributions  are 


I(x,y).  Then,  the  outputed  electric  current  on  each  detector  is 
[3] 


\\ 1  (*'y)dxdy 

!  At 


(*-1,2,  3,4) 


(1) 
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Fig . 1  Quadrant  detector  system 


Here, v  is  incident  light  frequency.  n  is  quantum  yield.  h 
isPlanck's  constant,  q  is  electron  charge.  Ak  is  the  faculae 
Qj-gaon  the  kth  detector .  In  this  way ,  the  transmission  signal  in 
thex  direction  is 


Sx=h  +  U~h~h 


(2) 


The  transmission  signal  in  the  y  direction  is 


Sy  —  +  h  h  *4 


(3) 


Angular  target  miss  amounts  are 


6,=  K-f  (4) 

*,-*-§*  (5) 


Here  S=il  +  i2  +  i3  +  i4.  K  is  a  proportionality  constant.  It 
is  obtained  through  zero  noise  signal  situations  and  represents 
tracking  system  angular  sensitivities. 

Besides  this,  angular  target  miss  amounts  can  also  be 
represented  as 


(6) 


J  2  6  V  2k  S/N 


V  2  stands  for  angular  target  miss  amounts  that  are  equal 
in  x  and  y  directions.  S/N  stands  for  signal  to  noise  ratio. 

This  clearly  shows  that  it  is  possible  to  reduce  tracking  error 
through  increasing  tracking  system  signal  to  noise  ratios. 
However,  with  regard  to  the  influences  of  atmospheric  turbulence, 
raising  signal  to  noise  ratior  is  to  no  avail. 


2  ATMOSPHERIC  TURBULENCE  NOISE  FORMS 


We  now  consider  the  noise  produced  by  turbulence  when  lasers 
are  propogated  in  the  atmosphere.  In  reality,  atmospheric 
turbulence  most  certainly  is  not  unifooa.  It  includes  a  good 
number  of  types  of  nonuniform  lens  media  [ 5 ] .  The  indices  of 
refraction  are  random  variables  associated  with  time. 

Fluctuations  in  indices  of  refraction  cause  light  waves  to 
produce  amplitude  and  phase  fluctuations  during  their 
propogation.  Under  certain  conditions,  these  fluctuations  are 
enough  to  cause  relatively  large  tracking  errors. 

Assume  the  angular  error  created  by  atmospheric  turbulence 
is  0t  and  that  angular  errors  led  to  by  phase  fluctuations  and 
light  intensity  fluctuations  are,  respectively,  0s  and  01. 

These  two  partial  errors  can  be  taken  to  be  statistically 
independent.  As  a  result,  we  have 


e  =  {e)+e)Vn 


(7) 


Below,  the  two  types  of  errors  are  discussed  respectively: 

(1)  Phase  Fluctuation  Due  to  phase  fluctuations  at 
reception  apertures,  beams  will  produce  expansion  and  vibration. 

vibration  and  reception  aperture  focal  length  then  determine 
instantaneous  beam  arrival  angle  variances  [6] 

2^ 


(a2)  = 


(8) 


DAp) 

k2p 1 


Here,  k=2rc/A.  A  is  wave  length.  p  is  reception  aperture 
diameter.  Ds(p)  is  phase  structure  function.  Option  is  made 
for  the  use  of  Von  Karman  turbulence  spectra.  With  regard  to 
spherical  waves,  it  is  possible  to  represent  them  as 


D,  (p)  =  1.091  C2fc2Lrp3/3[l -0.8048  (-£-  )'3  +  0.136  ( -^  )2]  (9) 


In  this  way,  arrival  angle  variances  can  be  written  as 

<aJ>  =  1.048C2(A^)-,,<Lr[l-0.8048  (o-^-  )>'3  +  0.136  (^-  )2]  (10) 
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Here,  c 2  is  atmospheric  turbulence  structure  constant.  Aa 

is  reception  aperture  area.  LT  is  target  range.  LO  is 
turbulence  outer  scale.  Establishment  of  the  equations  above 
assumes  reception  apertures  larger  than  4~XL  •  Root  mean 

sguare  electrical  currents  on  detectors  are  -T. 


h  V  4<  a2) 


(11) 


Signal  to  noise  ratios  are: 


S3 


s_ 

N 


(12) 


Therefore, 


VTe,  =  V  2<a2>  =  { 2.096  C*  (AX  )  ',/< 


Lt[  1-0.8048 


(_£_  ) «/»  +  0.136 
L0 


)J]  },/J 


(13) 


(2)  Light  Intensity  Fluctuations  If  we  assume  wave 
propogation  in  direction  z,  wave  fields  are  written  to  be 


(p  =  \J/exp  (-  ikz ) 


(14) 


Then,  i|i  can  be  represented  by  the  use  of  parabolic  form  equation 
[7] 


2ik3 t\j/  +  if/  +  2k1nxij/  =  0 


(15) 


Here  nl  =  n  -  <n>  is  index  of  refraction  undulation.  Ai  + 
Taking  extension  iriedia  and.  dividing  them  up  into  a  goo 
number  of  thin  layers  with  thickness  6z,  each  thin  layer  can  act 
as  a  random  phase  screen.  Each  screen  only  introduces  an 
independent  amount  of  phase  change  with  regard  to  waves 
Dropoqated  in  it,  and  amplitudes  do  not  change.  Amplitude 
changes  depend  on  a  good  number  of  6 z  diffraction  accumulations. 
As  far  as  this  time  is  concerned,  phase  correlation  functions  can 
be  represented  as 


(Pi  >  Pi ) 


<n,  (z,  p,)n,  (z  p2)  >dzdz' -k25,An  (p,-p2) 


(16) 


A4 


Here  p  is  the  two  dimensional  position  vector  (x,y).  An(p)  is 
the  correlation  function  associated  with  index  of  refraction 
fluctuations  on  plane  x-y.  It  is  also  capable  of  using  turbulent 
flow  spectrum  $n(k)  for  representation  as 


(p)  =2ttJ:2<51  jj  ( k i  k,  =0)exp  (f*  •  p)c t  k 


(17) 


As  to  processes  associated  with  flow  fields  from  one  screen  to 
another  screen,  it  is  possible  to  obtain  them  through  solving 
equation  (14).  Assuming  that  nl  =  0,  in  the  frequency  spectrum 
realm,  7(z,k)  is  a  two  dimensional  Fourier  transform  associated 
with  ^r(x,y,z).  In  that  cise,  equation  (14) 's  transformation  is 
made  into 

d.'F  (2.  Jt)--i  ~  ^(2,k)  (18) 

It  has  the  solution  . 

¥(*'  K)»y  u.»c)exp[-  lK-Yk~-  (19)  : 

In  this  way,  the  process  with  regard  to  wave  propogation  through 
6 z  =  z 1 -z  in  a  random  medium  is 

(a)  if>  (-,  p)==  ,  p)exp[  -ifc  n,  (p,  z  ")dz  "] ; 

Jo 

(b)  discrete  Fourier  transform  of  i|i  (z,p)  to  Y  (z,k); 

(C)  T  (z  ,k)  =  ^  (z,  K)exp[ — —  —  —  ]; 


(d)  inverse  discrete  Fourier  transform  of  7  (z‘-,k)  to  i|>  (z’~  ,p) 


Here — as  far  as  using  (z‘-,p)  to  represent  wave  field 
propogation  to  fields  associated  with  point  z  in  screens  is 
concerned — this  process  is  repeated  with  respect  to  all  phase 
screens  right  up  to  wave  fields  arriving  at  observation  planes. 
Finally,  the  squares  are  then  selected— that  is,  representing 
light  field  distributions  on  detectors.  Fig. 2  and  Fig. 3  are  then 
nothing  else  than  optical  intensity  distribution  graphs 
calculated  out  making  use  of  phase  screen  methods.  We  only 
select  data  from  middle  sections.  From  the  Fig.'s,  it  can  be 
seen  that  they  are  remote  Gaussian  distributions.  Making  /546 
use  of  this  light  intensity  distribution,  it  is  possible  to 
calculate  out  0t  .  Again,  making  use  of  equation  (7)  or  equation 
(13)  to  subtract  the  influences  of  phase  fluctuations,  that  is, 
it  is  possible  to  calculate  the  influences  of  light  intensity 
fluctuations  on  angular  precisions. 


Fig. 2  Intensity  distribution 


Fig. 3  Intensity  distribution 


3  CALCULATION  RESULTS  AND  CONCLUSIONS 

We  considered  cases  of  horizontal  laser  propogation  through 
the  atmosphere.  The  transmiting  laser  wave  length  is  1.06Mm. 
Divergence  angles  are  1.5mrad.  Making  use  of  Von  Kantian 
turbulence  spectra,  turbulence  outside  scales  Lo  are  2m.  Inside 
scales  are  0.01m.  Reception  apertures  are  0.4m.  As  far  as  the 
influences  of  phase  fluctuations  on  laser  tracking  systems  are 
concerned,  it  is  possible  to  make  use  of  equation  (12)  to  carry 
out  calculations.  As  to  the  influences  of  intensity 
fluctuations,  we  assume  that  tracking  systems  have  already  been 
aligned  on  targets.  Phase  screen  methods  are  used  to  calculate 
light  intensity  distributions  on  detectors.  During  calculations, 
we  make  use  of  512x512  screens.  Screen  intervals  are  100m. 

First  of  all,  calculations  are  made  of  beam  field  distributions 


at? 


through  the  first  screen.  On  the  basis  of  already  known  tracking 
distances,  precise  determinations  are  made  for  the  number  of 
screens  through  which  light  fields  pass  right  up  to  obervation 
planes.  Finally,  calculations  are  made  of  light  intensity 
distributions  on  observation  planes.  Giving  consideration  to 
limitations  on  reception  apertures,  we  only  select  light  field 
distribution  data  inside  apertures.  Making  use  of  equations  (4) 
and  (5),  calculations  are  done  of  the  influences  of  atmospheric 
turbulence  on  laser  tracking  angular  precisions.  Here,  the 
proportionality  constant  K  =  1.473xl0-4/f  rad.  Then,  use  is  made 
of  equation  (7)  to  calculate  out  the  influences  of  light 
intensity  fluctuations.  In  calculations,  first  of  all,  it  is 
assumed  that  transmission  distances  are  fixed  (distance  is  5km). 
Turbulence  intensities  change.  Calculation  results  are  shown  in 
Fig. 4.  After  that,  it  is  assumed  that  turbulent  flow  intensities 
are  fixed.  Transmission  distances  are  varied.  Results  are  seen 
in  Fig. 5.  Table  1  sets  out  calculation  results  for  i0,  and  0j 
under  these  conditions.  From  Table  1,  it  is  possible  to  know 
that  tracking  errors  created  by  intensity  fluctuations  play  the 

Table  1  Tracking  error  induced  by  atmospheric  turbulence 


CM 

1  x  10'u 

IX  10 

l-M 

LT  /km 

0(//*rad 

0,/  p  rad 

OJfirad 

$,/ fiT&d 

1 

1.25 

4.62 

3.96 

11.22 

2 

1.77 

6.03 

5.60 

13.75 

3 

2.17 

7.39 

6.86 

15.92 

4 

2.50 

8.67 

7.92 

17.75 

5 

2.80 

9.90 

8.85 

19.21 

6 

3.07 

11.06 

9.70 

20.33 

7 

3.31 

12.15 

10.48 

21.10 

8 

3.54 

13.18 

11.20 

21.51 

9 

3.77 

14.15 

11.89 

21.57 

10 

3.96 

15.05 

12.52 

21.27 

/547 

dominant  role.  Tracking  errors  led  to  by  phase  fluctuations  are- 
— when  turbulent  flows  are  relatively  weak  approximately  one 
third  of  intensity  fluctuations.  Moreover,  when  turbulent  flow 
strengths  are  medium,  the  difference  between  the  two  is  somewhat 
reduced.  The  farther  tracking  distances  are,  the  smaller  the 
difference  is.  This  phenomenon  may  possibly  be  created  by 
saturation  effects  [6]  associated  with  intensity  fluctuations. 
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Fig. 4  Tracking  error  vs  turbulence 


Fig. 5  Tracking  error  vs  range 


On  the  basis  of  our  observations ,  as  far  as  turbulent  flow 
strengths  near  the  earth's  surface  are  concerned,  they  are  weak 
when  under  10-16m-2/3.  During  normal  clear  weather,  they  are 
around  10-14m-2/3.  When  they  are  strong,  they  have  the 
capability  of  reaching  above  10-13m-2/3.  From  these  Fig.'s,  it 
is  possible  to  see  that — in  cases  where  turbulent  flows  are  weak- 
-for  distances  inside  10km,  tracking  errors  given  rise  to  by 
atmospheric  turbulence  do  not  exceed  5/irad.  In  situations  where 
turbulence  is  moderate,  tracking  errors  reach  25/xrad.  If 
turbulence  is  even  stronger,  then,  at  close  ranges,  they  are  also 


capable  of  exceeding  25jurad. 

Speaking  in  terms  of  aerial  target  tracking,  the  higher 
altitudes  are,  the  weaker  turbulence  is,  and  the  smaller  tracking 
errors  are.  However,  the  closer  it  is  to  the  earth’s  surface, 
the  stronger  turbulence  is.  This  is  the  case  for  no  other  reason 
than  that — when  angles  of  elevation  are  low — it  is  difficult  to 
track.  As  to  tracking  errors  which  are  created  by  turbulence  in 
this  case,  we  will  discuss  them  again  later. 

Acknowledgements:  Our  thanks  to  Comrade  Wang  Yingjian  for 
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EFFECTS  OF  COATING  THICKNESS  NON- UNIFORMITY  ON 
FAR-FIELD  PERFORMANCE  OF  BEAM  INDICATOR 
Xiong  Shengming  Zhang  Yundong  Wang  Yongzhong 

ABSTRACT  Optical  thin  films  in  DF  (deuterium  fluoride)  high 
power  laser  system  beam  directors  are  discussed.  Use  is  made  of 
polarization  ray  trace  to  analyze  polarization  state  changes  in 
main  lasers  and  beacons  at  director  pupils.  Calculations  are 
made  of  the  influences  of  coating  thickness  nonuniformities  on 
far  field  performance. 


KEY  WORDS  High  power  laser  Beam  director  Multilayer 
coastings  Coating  thickness  nonuniformity  Far  field  performance 


30 


0  INTRODUCTION 


High  power  laser  systems  require  that — in  optical 
components — use  be  made  of  muliple  layer  optical  coatings  to  make 
their  high  efficiency  output  laser  energies.  On  optical  surfaces 
which  have  curvatures,  incident  angles  across  optical  apertures 
continuously  change.  In  particular,  within  incident  angle  ranges 
associated  with  beam  expander  device  reflector  surfaces, 
reflectance  changes  (Rs  and  Rp)  are  not  obvious.  However,  as 
far  as  across  apertures  on  the  same  reflectors  is  concerned, 
incidence  angles  continuously  change  and  optical  coating 
reflectances  as  well  as  reflection  phases  change  following  along 
with  wave  length.  On  entire  mirror  surfaces,  phase  extensions  A 
(0s  -  0p)  change  following  along  with  incidence  angles  and  wave 
lengths.  Small  coating  thickness  changes  will  still  lead  to 
obvious  A  changes  [1,2]. 

In  large  aperture  laser  emission  systems,  as  far  as  sensing 
devices  which  are  sensitive  to  polarization  are  concerned — for 
example,  interferometers,  wavefront  sensors,  tracking  devices, 
and  so  on — in  measurements  of  such  data  with  regard  to 
wavefronts,  it  is  necessary  with  respect  to  entire  light  path 
systems  to  complete  polarizaton  analyses  of  them  from  beginning 
to  end.  As  far  as  light  paths  associated  with  multiple  layer 
optical  coating  optical  components  and  possessing  coating 
thickness  changes  are  concerned,  coating  thickness  changes  are 
capable  of  giving  rise  to  wavefront  changes  associated  with  peak 
and  valley  differences  which  reach  a  wave  length.  If,  in  regard 
to  wavefront  detection,  use  is  made  of  light  wave  and  emitted 
laser  light  wave  lengths  which  are  not  the  same,  as  a  result, 
wavefronts  measured  by  wavefront  sensor  devices  are  not  the  same 
in  wavefront  errors  [3]  associated  with  emitted  laser  light 
undergoing  the  same  light  path.  Detected  wavefronts  and  emitted 
laser  light  wavefronts  change  following  along  with  incidence 
angle  and  laser  operating  time  developments. 

This  article  discusses  optical  coatings  used  by  DF  high 
power  laser  systems.  It  uses  polarized  light  ray  trace  analysis 
on  pupil  polarization  changes  and  calculates  the  influences  of 
coating  thickness  nonuniformities  on  far  field  performance. 


1  REFLECTOR  FILM  SYSTEMS  AND  CHARACTERISITCS 


Based  on  overall  optical  system  requirements,  we  opt  for  the 
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use  of  metallic  media  strengthened  film  systems.  The  design  of 
the  entire  film  system  is  Sub/CrAg  (PbF2/ZnSe) 3/Air .  In  it,/593 

Ag  thickness  is  lOOnm.  PbF2  and  ZnSe  optical  thicknesses  are  1/4 
of  3.8/im.  Fig.l  gives  visible  and  infrared  theoretical  spectrum 
curves  associated  with  this  type  of  coating  system  design.  On 
the  basis  of  currently  available  film  plating  equipment  as  well 
as  film  plating  technology,  in  beam  directors,  $600  main 
reflector  coating  thickness  change  is  5%-6%.  Secondary  reflector 
film  thickness  change  is  adopted  as  3%.  Coating  thickness 
distributions  have  their  source  in  the  calculations  of  Reference 
[2].  Fig. 2  is  beam  director  main  reflector  coating  thickness 
distribution  forms.  In  high  power  laser  transmission  systems, 
other  optical  component  coat'ng  thickness  distribution  forms  are 
analogous  to  this. 


Fig.l  Theoretical  reflectance  spectral  for  coating  Sub/CrAg 
( PbF2 /ZnSe ) 3 /Air 


3-a- 


Fig. 2  Coating  thickness  variation  model  of  primary  mirror  of 
director 
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Fig. 3  Phase  variation  due  to  nonuniformity  for  standard 
quarterwave  design 


In  multiple  layer  coating  systems,  small  coating  thickness 
variations  (CTV)  will  give  rise  to  large  equivalent  surface  form 
errors.  Surface  form  errors  which  are  given  rise  to  by  coating 
thickness  variations  are  produced  by  reflected  light  phases 
because  reflected  light  phases  are  functions  of  coating 
thickness.  Fig. 3  represents  reflected  light  phase  changes  given 
rise  to  by  coating  thickness  variations.  In  the  Fig.,  angle  of 
incidence  is  15°.  Reflectance  of  3.8Aim  wave  length  laser  light 
is  99.9%.  Below,  we  will  investigate  the  sensitivities  of  system 
properties  to  these  changes. 

In  studies,  it  is  assumed  that  coating  thickness  varies 
along  with  radial  direction  changes  and  rotations  are 
symmetrical.  Coating  thickness  variations  are  dependent  on 
optical  component  size,  the  geometrical  dimensions  of  coating 
plating  machines,  and  other  plating  techniques.  In  situations 
where  multiple  high  reflectance  coatings  have  coating  thickness 
variations,  reflectance  changes  are  very  small.  However, 
reflected  light  phases  are  capable  of  changing  several  hundred 
degrees . 


2  POLARI Z AT I ON  ANALYS I S 


Before  we  do  some  work,  there  is  still  no  polarization 
analysis  software  inside  China  which  is  capable  of  handling  thin 
films  that  possess  coating  thickness  variations  associated  with 


symmetrical  rotations .  As  far  as  our  writing  of  polarized  light 
ray  trace  computer  software  is  concerned,  use  was  made  of 
analytical  methods  associated  with  Jones  calculations  [4]  and 
Eugene  Waluschk's  discussions  [5].  This  software  studies  the 
taking  of  each  polarized  component  (each  component  in  a  light 
path)  and  dividing  it  into  a  combination  of  parts  that  are 
polarization  devices  and  delay  devices.  Polarization  device 
light  axses  are  parallel  to  the  beam  and  the  plane  s  of  the 
optical  surface  affecting  it.  Delay  device  light  axses  are 
parallel  to  the  plane  p. 

Reflector  polarization  device  matrix  primary  transmission 
rates  are  Rs  =  r),  and  Rp  =  r*  .  Rs  and  Rp  are, 

respectively,  s-  an!  p-  component  reflectances.  <ps  and  0p  are 
plated  coating  reflector  s-  and  p-  component  reflection  phases. 
Polarization  device  matrices  are  complex  numbers.  Generally,  rpi 

*  rsi  .  Selecting  the  incident  surface  included  angle  /594 
oetween  reflectors  to  be  0i,  reflection  matrices  associated  with 
various  reflectors  are 


0 

r,  e'9' 


Conversion  matrices  associated  with  light  rays  entering  into 
reflector  corrdinates  are 


sinfl, 

COS0j 
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Reflector  coordinate  reconversion  matrices  are 


Jones  matrix  calculations  obtain  an  array  of  vectors.  It  is 
composed  of  two  polarization  states  which  are  mutually 
perpendicular  to  each  other  with  amplitude  and  phase  being, 
respectively.  Ax,  8,45 
x.  Ay,  <py 


In  order  to  precisely  determine  the  influences  of  thin  films 
on  surfaces  with  curvatures  as  well  as  the  influences  of  coating 
thickness  variations  on  reflector  accumulations,  it  is  necessary 
to  calculate  Jones  vectors  associated  with  different  aperture 
positions  across  laser  systems.  We  are  able  to  recognize  that 
laser  system  pupils  are  two  types  of  different  apertures.  One 
type  is  S  polarized  light  apertures.  The  other  type  is  P 
polarized  apertures.  Light  wave  polarization  states  inside  the 
two  apertures  are  mutually  perpendicular.  As  far  as  carrying  out 
calculations  with  regard  to  variations  in  amplitude  and  phase 
associated  with  the  two  apertures  is  concerned,  two  dimensional 
Fourier  transforms  are  then  carried  out  on  each  aperture's  wave 
surface.  Far  field  strength  distributions  are  then  obtained. 

The  two  polarization  state  far  field  strength  distributions  are 
then  added  together.  Eliminating  2,  overall  strength 
distributions  are  then  obtained.  The  ratios  between  far  field 


strength  distribution  center  faculae  intensities  and  ideal  system 
central  faculae  intensities  are  Strehl  ratios. 


3  CALCULATION  RESULTS 


On  reflectors,  multiple  layer  optical  coatings  will  give 
rise  to  reflected  light  phase  delay  variations,  thereby  leading 
to  obvious  polarization  and  wave  front  changes.  For  this  reason, 
we  calculated  polarization  state  variations  and  far  field 
strength  distributions  in  pupils  associated  with  different  wave  r 

lengths  of  beam  director  light  waves.  In  calculations,  main 
reflector  coating  thickness  variations  are  selected  as  5%. 

Secondary  mirror  coating  thickness  variations  are  selected  as  3%. 

Coating  systems  are  the  same  as  strengthened  metallic  media 
reflection  coatings.  In  calculations,  consideration  has  not  yet 
been  given  to  chromatic  dispersion  with  regard  to  different  wave 
lengths  and  media  coating  layers.  Fig. 4  represents  horizontal 

ray  polarized  incident  3.8/um  wave  length  laser  light  in  two  cases 
of  pupil  polarization  state  changes — beam  director  reflectors 
with  and  without  coating  thickness  variations.  Fig. 5  is  pupil 
polarization  state  changes  associated  with  0.6328/im  beacon  light. 

Fig. 6  and  Fig. 7  are,  respectively,  far  field  intensity 
distributions  at  beam  directors  associated  with  3.8/xm  intensity 
light  and  0.6328/im  becon  light.  Table  1  summarizes 


Fig. 4  The  polarization  state  of  3.8/xm  laser  at  the  pupil  of 
director  (a)  Without  coating  thickness  variation  (b)  Coating 
thickness  variation  (main  mirror  5%,  secondary  mirror  3%)  / 595 


the  influences  of  coating  thickness  variations  on  three  main 
spectral  lines  associated  with  DF  laser  light  transmission  and 
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beacon  wave  length  light  wave  far  field  Strehl  ratios . 


(a)  (b) 


Fig5 .  The  polarization  state  of  0.6328/im  laser  at  the  pupil  of 
director  (a)  Without  coating  thickness  variation  (b)  Coating 
thickness  variation  (main  mirror  5%,  secondary  mirror  3%) 


Fig. 6  3.8/xm  laser  far-field  intensity  distribution 


3? 


///« 


Fig. 7  0.6328/im  laser  far-field  intensity  distribution 
Table  1  Far  field  performance  effects  due  to  coating  thickness  variation  (CTV) 


wavelength  /  Atm 

0.6328 

3.80 

3.716 

4.041 

no  CTV 

0.9541 

0.9941 

0.9941 

0.9941 

5%  CTV  of  primary  mirror 

3%  CTV  of  secondary  mirror 

0.8504 

0.9607 

0.9498 

.  0.9783 

4  CONCLUSIONS 

if  laser  systems  do  not  carry  out  overall  wavefront 
calibration,  polarization  effects  given  rise  to  by  thin  coatings 
S?i l  havi  obvious  influences  on  wavefronts  and  far  field 
erfoSce  fc  io  multiple  layer  coating  interference  effects, 
wavefront  errors  undergone  by  main  laser  and  beacon  light  are  not 
the  Se!  booking  at  the  calculation  results  above,  influences 
of  coating  thickness  variations  on  beacon  light  polarization  and 
far  field  performance  are  much  more  severe  than  high  power 
lasers61  Using  the  calculations  above,  consideration  was  only 
made  of  two  director  reflectors.  In  the  calculations  above, 
option  was  made  for  the  use  of  coating  system  polarizations  that 
are  relatively  small.  0.6328Mm  wave  length  beacon  light  far 

ratios  then  fall  to  0.85.  This  already  approaches 
the“raSScnaf Rayle^h  criterion.  In  situations  equivalent  to 
5%  coating  thickness  variation,  reflectors  are  not  able  to  /596 

eXiStIf°conl!d:r™is  given  to  all  reflectors  associated  with 
nrvHral  svstems — in  particular,  those  refraction  transfer 
reflectors  where  coating  thickness  variations  also  exist-- it  will 
ltld  to  severe  polarization  effects,  producing  image  distortion 
i  oo  n-F  fnrus  ril  Far  field  performance  will  drop  a  step 
ISrtie?  Lawsuit;  thSre  is  a  need  to:  (1)  optimise  coating 
system1" design,  control  phase  delay  acr< ass 

coating  thickness  variations,  as  well  as.the. ontimize 
inriHpnt  angles  on  polarization  and  wavefronts,  (2)  optimiz 
cSrting  plating  techniques,  reducing  coating  thickness  errors. 
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(3)  go  through  entire  system  performance  calculations  to 
precisely  determine  coating  thickness  errors;  (4)  match  up 

residual  wavefront  errors  associated  with  high  power  lasers  and 
Deacon  light. 
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